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Abstract 

We study spinor field theories as an origin to induce space-time evo- 
lution. Self-interacting spinor fields with canonical and non-canonical ki- 
netic terms are considered in a Friedman- Robertson- Walker universe. The 
deceleration parameter is calculated by solving the equation of motion and 
the Friedman equation, simultaneously. It is shown that the spinor fields 
can accelerate and decelerate the universe expansion. To construct real- 
istic models we discuss the contributions from the dynamical symmetry 
breaking. 

1 Introduction 

The acceleration periods of the universe evolution are one of the most important 
problems to be explored in astroparticle physics today. Though the cosmological 
constant gives a simple solution to accelerate the universe expansion, a time 
dependent source is necessary to stop the inflation at the early universe and 
restart it at the present low energy scale. There is no simple candidate to derive 
the accelerated expansion of the universe in the Standard Model contents for 
particle physics. The problem calls for considerations to new field theory models 
which may describe the acceleration periods of the cosmic expansion. 

Various models beyond the Standard Model have been studied to include 
candidates which accelerate the expansion of the universe. The simplest model 
can be constructed by an inflaton like scalar field which is easy to treat in a 
curved spacetime. A non-vanishing potential energy for the scalar field can 
induce the exponential expansion of the universe. On the other hand, a spinor 
field also has played an important role as a gravitational source. An ordinary 
free fermion behaves as a matter field and decelerates the universe expansion. 
Much interest has been demonstrated in order to extend the model for the spinor 
field in recent years. 

The study of spinor fields in curved spacetime has a long history. The Dirac 
equation was investigated for massless spinor fields in curved space-time more 
than 50 years ago [T]. The Einstein-Dirac equations were solved for a massive 
spinor field in a special anisotropic spacetime, Bianchi type I universe with a 
cosmological constant [2J. The possibility to induce a primordial inflation [3] 
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and the current expansion in a self-interaction spinor field [J] has recently been 
pointed out. Similar self-interacting fields have been discussed as a kind of 
Inflaton [5] and Quintom [5J [7] . They have been studied in a Bianchi type I 
framework [5J [SJ HH] and other anisotropic space-time [TTJ [T2] . Non-standard 
spinors are considered as a source for the current expansion in Refs. [T3l [14j 
[T51 H HO [H HI H3 HJ. There is a possibility to cause late-time cosmic 
acceleration by a spinor field through a coupling with the Brans-Dicke scalar 
field [221 123] , a form invariance transformation |24| and a non-minimal curvature 
coupling [2"5] . 

The class of dark energy models has been proposed in a scalar field with 
non-canonical kinetic terms, such as models named k-inflation [26l [27] or k- 
essence [551 123 ISO]- A similar extension is possible for a spinor field. The 
spinor counterparts of k-essence for scalar field are called f-essence [3TJ [32] and 
g-essence [33]. We would like to continue the study of such models on self- 
interacting spinor fields with non-canonical kinetic terms as well. It should be 
noted that models with the square of an usual Dirac Lagrangian may be special 
cases of f-essence scenario [3H [35] . 

The scalar invariant constructed from two spinor fields dynamically develops 
a non- vanishing value in QCD like theories [351 [37]. Then the chiral symmetry 
for the spinor fields is broken. It has played an essential role for the evolution 
of the universe. Only a little work has been done to evaluate the dynamical 
symmetry breaking in a non-static spacetime. In de-Sitter space it is possible 
to calculate the expectation value for the scalar invariant based on the maximal 
symmetry of the space [35J [35]. However, we would like to find solution with 
early- and late-time acceleration periods. For this purpose it is not avoidable 
to develop different procedure. 

In this paper a spinor field is assumed as a gravitational source. We confine 
ourselves, for simplicity, in a Friedman-Robertson- Walker (FRW) metric and 
evaluate the corresponding space-time evolution. In Sec. 2 we review the space- 
time evolution induced by a spinor field with an ordinary canonical kinetic term. 
A self-interacting spinor field is considered in curved space-time. Evaluating the 
equation of motion for the spinor field and Einstein's field equations, we show 
typical behaviors of the Hubble and deceleration parameters. It is quite general 
to introduce a non-canonical kinetic term in an effective model. We assume a 
simple form for the non-canonical kinetic term and evaluate the influences on 
the space-time evolution in Sec. 3. In Sec. 4 we set the critical scale where 
the dynamical symmetry breaking takes place by hand and discuss the critical 
behavior caused by the dynamical symmetry breaking. The space-time evolu- 
tion may be suddenly modified by a first order phase transition. Finally, some 
concluding remarks are given in Sec. 5. 

2 Spinor fields with a canonical kinetic term 

Before we study special features of the non-canonical kinetic term for a spinor 
field, it is more instructive to discuss the space-time evolution induced by a 
spinor field with an ordinary canonical kinetic term. The space-time evolution 
is evaluated by solving the equation of motion for the spinor field and Einstein's 
field equations, simultaneously. The variations of the action give rise to these 
equations. 
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Here we start from the action, 



S = J d^x4~g {L g + C A + C D ) , (1) 



where C g is the Einstein-Hilbert Lagrangian and Ca describes a cosmological 
constant term 

16^G R,Ca = ~8^G" 
with R the curvature scalar and G the Newtonian constant. The Dirac La- 
grangian Cd is given by 



A, = — ^a = -^ A > ( 2 ) 



C d = \ W" D ^ ~ {D^)T^} - miM> - V(^), (3) 
where the covariant derivatives for spinor fields, ip and ip, are defined by 

= drf - (4) 

and 

D$ = drf - ^0 M . (5) 
The spinor connection, is written by the vierbein, e M a , 

% = --gpaF"^ - e^e^PT 4 , (6) 

and the Dirac matrices, r M , are generalized in a curved spacetime, 

P 1 = e" a7 °. (7) 

Here we suppose that the main contribution as a gravitational source comes 
from a non- vanishing expectation value for a scalar invariant, ipip. The potential 
V(ipip) can be written as a function of a composite operator ipip. We noticed that 
only scalar invariants can develop non- vanishing vacuum expectation values to 
keep the Lorentz structures in Minkowski spacetime. It is not valid in a curved 
space-time. Thus the above assumption should be modified in a strongly curved 
space-time. 

The equations of motion for ip and ip are given by the variations of the action 
(JT|) with respect to ip and ip, respectively, 



dV 

dip 

and 



i(D^)T» - mil - — = 0, (8) 



iT^Daip - mtp - 4t = 0. (9) 
dip 

The variation of the action ([1]) with respect to the vierbein gives the Einstein's 
field equations, 

1 

2" 

where the energy-momentum tensor, T^, is given by 



- ■zRg^ + A 5m „ = -8vrG(7V), (10) 



V = l - [ipT^DviP + i>T v D^iP - {D y {p)T^ - (D M V0r„V] 

-g^C D . (11) 
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The space-time evolution is found by solving Eqs. flSJ), © and ([TU| with the 
help of the above energy-momentum form. 

We consider the case that the spinor field has a homogeneous and isotropic 
distribution and look for the solution in a FRW universe with a flatly spatial 
part. The FRW universe is defined by the metric 

ds 2 = di 2 -a 2 {t)(dx 2 + dy 2 + dz 2 ). (12) 

The space-time evolution is described by the time development of the scale 
factor a(t). Under the FRW metric, (fT2")) . the generalized Dirac matrices are 
given by 

r° = 7°, ^ = -^7*, (13) 
a(t) 

and the spinor connection reads 

r> = 0, a = id(i) 7 l 7°. (14) 

Since the spinor field is assumed to be homogeneous and isotropic, it is 
written as only a function of time. In this case Eqs. (JHJ and (JSJ) are simplified 
to be 

d tp + ^Htp - imtPj - i— 7° = 0, (15) 
2 dtp 

and 

doip + ^HiP + rfmip + i^ ^ = 0, (16) 
2 dtp 

where H is the Hubble parameter, H = a(t)/a(t). The derivative of the poten- 
tial, V(tptp), can be written as 

From Eqs. (TTSl) . (fTc?)) and (fTTf the potential dependence can be eliminated. 
Therefore the composite operator, ipip, satisfies a simple equation of motion, 

(Sb + aff)^v = o. (18) 



This equation has a trivial solution, ipip — 0. A non-trivial solution of Eq. (|18p 
is 

a 3 (<) 

where C is an arbitrary constant parameter fixed by an initial condition [3J 01 E] ■ 
The composite operator, ipip, develops as an ordinary matter field. We noticed 
that the solution (TH))) is obtained at the classical limit. It may be modified 
by radiative corrections for the potential and the energy-momentum tensor. 
It should be also noticed that the solution is found for a spatially constant, 
ipip. However, we consider that the solution has something essential for the 
development of the composite operator and discuss the evolution of the space- 
time under the solution (fl^|. 

In the FRW universe the energy-momentum tensor, T{f, has diagonal form. 
It can be parametrized as 

= diag(p, —p, —p, —p), (20) 
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where p corresponds to energy density and p pressure. Thus the Einstein's field 
equations ([TUl) read 

'a(ty " 



a(t) 



A = 8irGp, (21) 



and 

l)-(f) 3 + A — <»> 

From Eq. (jTTJ) with the equation of motion (fT5|) and (fTo]) the energy density, 
p, and the pressure, p, can be expressed as functions of the scalar invariant, ipip, 

p = m(^) + (V(^)}, (23) 

and 

P=-(V(^j)) + (V'(^)-^}. (24) 

An explicit form of the potential is necessary to find the typical behavior of the 
space-time. It is assumed that the potential can be expanded in terms of the 
scalar invariant, 

oo 

(y(^)) = ^a„(^) 2 ", (25) 

n=l 

where a n corresponds to a coupling constant for multi-fermion interactions. To 
avoid an unexpected divergence we consider only a positive n in Eq. (|2"5)l . Since 
the scalar invariant, ipi/j, is a dimension 3 operator, the potential (|23|) is not 
rcnormalizable. We regard the potential as a low energy effective one comming 
from a more fundamental theory at a high energy scale. It is expected that 
the higher order terms of ip^ are suppressed at low energy. Below, we cut 
the summation at n max . We noted that the model is a kind of four-fermion 
interaction models for n max = 1 which is often used as a low energy effective 
model of QCD [31137]. 

Hence, the energy density, p, and the pressure, p, read 

p = m{W>)+ £ «n(#) 2n , (26) 

and 

"max 

p= a n {2n-l){^) 2n . (27) 

n=l 

For a trivial solution of Eq. (p~8|) . (ipip) = 0, the right hand sides of the Einstein's 
field equations, (|2"Tj) and (|2"2"|) vanish, which simply describes a well known de 
Sitte universe. We can recognize that the solution corresponds to lepton and 
heavy quark fields. In our universe finite vacuum expectation values are not 
observed for the lepton and the heavy quark bilinears. These spinor fields have 
no role for the evolution of the universe. However, it is known that the composite 
operator, ipip, develops non-vanishing value for up and down quarks in the 
current universe. Thus we evaluate the space-time evolution under the solution 
(Ti"T)l) . Inserting the solution (TT9|) into Eqs. (|2l)|) and (|27p . we rewrite the Einstein's 
field equations, (|2"Tj) and (|2"2"j) . as 



o(t) 



£j Uma x £j2n 
l_ \ ^ 



a 3 (t) ^ n a en (t) 



(28) 



5 



and 



a(t) 



8ttG 



d(t) 
o(t) 

"ma rc 

2 a„(2n-l) 



A 



G 2 



a 6 "(i) 



(29) 



As is known, these equations are not independent. Equation (1291) can be derived 
from Eq. (|2"5]l plus the Bianchi identity. These equations are called Friedman 
equations. 

The Hubble parameter ,H, then is given by 



H 2 = 



A 



IttG 



C 



n mQj; 

E 

n=l 



g 2 



-, 6n 



(*) 



(30) 



The cosmological constant and the mass terms mainly contribute to the right- 
hand side of Eq. (13U1) at the low energy limit, a(t) — > oo. As is well-known, the 
composite operator, tpip, develops a negative expectation value for up and down 
quarks in the current hadronic phase. It means that the parameter, G, has to 
be negative for the light quarks. Thus the mass term in Eq. ([50)1 is negative in 
the standard model of particle physics. On the other hand, the left-hand side 
of Eq. ([3"0f has possessed a non-negative value. The cosmological constant is 
necessary to take the low energy limit as, a(t) — > oo. Upper limit for the scale 
factor exists without the cosmological constant. It should be noted that the 
parameter G is assumed to be positive in Refs. [UE]. 
The deceleration parameter is written as 



a(t) 1 , 

where u shows the equation of state of the universe, 

-A + 8ttGp 



(31) 



A + 8irGp 
-A + 8?rG 



Ea n (2n — 1) _ , . 
y 'a 6n (t) 



A + 8ttG 



G 



E 



C 



2n 



a 6n (t) 



(32) 



An acceleration period is realized for a negative q. At the low energy limit, 
a(t) — > oo, it diverges for A = and approaches a negative unity for a non- 
vanishing A. 

Evaluating Eqs. (f3T)f and (|3"Tj) numerically, we show features of the spinor 
field as a gravitational source. In numerical calculations all the mass scales are 
normalized by the fermion mass to. In Fig. [T] typical behavior for the Hubble 
and the deceleration parameters are shown as a function of the scale factor a 
for a four- fermion interaction model without cosmological constant. An eight- 
fermion interaction is introduced in Fig. [5] The parameter, G, in Eq. (|30[) 
is negative for QCD-like theories. The equation (|3T)|) has a real solution for 
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Figure 1: Behavior of the Hubble and the deceleration parameters for A 
C = — 1, aim 2 = 5, 10, 15 and n max = 1. 
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Figure 2: Behavior of the Hubble and the deceleration parameters for 
A/(87rGm 4 ) = 0.1, C = -1, aim 2 = 5, 10, 15 and n max = 1. 



the Hubble parameter, only if the coupling constant for the higher dimensional 
operator, a nmax , takes a positive value. As is seen in Fig. Q](a), the right-hand 
side of Eq. (1501 develops a negative value for a larger a. Since the Friedman 
equation has no real solution, the parameter range for the dashed line is ruled 
out. There is an upper limit for the scale factor. In Fig. Q] (b) we observe that 
the deceleration parameter, q, blows up at the upper limit. 

In Fig. [5] we introduce a cosmological constant term for a four-fermion inter- 
action model. Since the cosmological constant term makes the right-hand side of 
Eq. (I5D1 shift, it is possible to cancel a negative contribution from the fermion 
self interaction terms. In other words, a cosmological constant like gravita- 
tional source is necessary to realize an open universe with a non- vanishing tptp. 
In Fig. [2] (b) the deceleration parameter converges to unity for a large a limit. 
Thus the late time universe can be accelerated by the cosmological constant. 

For a small a higher dimensional operators modify the behavior of the Hubble 
and the deceleration parameters. Interesting behavior is found if one of a n is 
negative. In Fig. [3] we draw the behavior of the Hubble parameter for a negative 
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Figure 3: Behavior of the Hubble parameter for A/(87tGto 4 ) = 0,0.1, C = — 1, 
aim 2 — 10, a2i)i 8 = —0.6, —0.8, —1, a^m 14 = 0.03 and n max = 3. 
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Figure 4: Behavior of the deceleration parameter for C 
q^to 8 = —0.6, —0.8, —1, c^m 14 = 0.03 and n max = 3. 
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eight-fermion coupling, a2 < 0. We also introduce a positive twelve-fermion 
coupling, Qi3 > to stabilize the potential. A negative ai modifies the behavior 
of Hubble parameter in Fig. [3] (a). As is shown in Fig. [3] (b) the square of the 
Hubble parameter becomes negative without cosmological constant. Thus the 
Friedman equation has no real solution on the parameter range for the dashed 
line. This again means that the universe has an upper limit. The deceleration 
parameter blows up at the scale in Fig. 2] (a). In the case a^jr? = — 1 we find 
an acceleration period, q < 0, without any cosmological constant. In Fig. 2] (b) 
it is found that the late time universe can also be accelerated by introducing 
cosmological constant. 



3 Non-canonical kinetic term 

In the previous section we only consider scalar type multi-fermion interactions. 
It is assumed that such interactions are introduced in a field theory as low energy 
effective model stemming from a more fundamental theory at high energy scales. 
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Varieties of interactions which keep the actual symmetry for the particle physics 
can be generated. A kinetic term for fermion may also be modified at low energy 
scale. Here we investigate a class of fermion field theories which are described 
by the Dirac Lagrangian with non-canonical kinetic terms, 

+m$i/>-V$ip), (33) 

where the function f(ipip) modifies the kinetic term. If we set f(ipip) = 1, the 
Lagrangian is reproduced. Here we consider a more general case, where the 
function f("ipip) is given by a function of only a scalar invariant -tfjip- 

The non-canonical kinetic terms modify the equations of motion for fermion 
fields, tp and -0, to 

-i/(D^)r" + mi> - V'rp = 0, (34) 



and 



+i.fT tl D^ + mip- V'tp = 0, (35) 



where we write 



d{ipip) 

Derivatives of the potential, V(ipip), and the function /('tpip) can be written as 
dV t/' 7 dV t/' / d -f t> J d f f> I (1-7\ 

M, =V d^ = V ■*> dj = f ¥ = / (3?) 

From Eqs. (fM)) and ([33)) with Eqs. (I3T1) we obtain a simple equation of motion 
for the composite operator, ipip, 

[(/ + #/l% + 3/F]# = 0, (38) 

There is a trivial solution, ipip = 0. To find a non-trivial solution we expand the 
function f(ijjip) in terms of the scalar invariant, 

= X^n^) 2 ". (39) 

Thus the non-trivial solution of Eq. is found to be 

/(^O&f = E MW) 2n+1 = (40) 

where C is a constant parameter. 

We numerically solve Eq. (|40|) for a fixed a(t = to), (i>ip)(t — t ) and (3 n . We 
draw typical behavior of ('tptp) in Fig.[5]with ( , ipip)(t = to) = m 3 and a(t = to) = 
mT l . The term, f3i(ipip) 2 , changes the slope for a small a. If the coefficient, /3i, 
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Figure 5: Behavior of the expectation value (ipip) f° r (V'V'H^ — ^o) — w 3 ; a(t 
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is negative, a positive fa is necessary to satisfy Eq. (|4"0")) for a large In 
this case a more gradual slope is observed for a larger fa. The term, j3 n {i)i>) 2n , 
with a negative n diverges at the limit, t/'V' ~~ * 0- Non-canonical kinetic term 

remains at a late time universe. Thus the term, /3 i (-0-0) 2 , is less interesting. 

In Fig.[5](c) a negative raises up a value, |^|, for a large a. 

In the FRW universe the Einstein's field equations are given by Eqs. (f2"Tj) and 
(|2"2"1) . For the Lagrangian (155)) the energy-momentum tensor, T^, is modified 
to be 

-(-D M V>)r„^] - g^Cr,, 

(41) 

From Eq. (|4"T|) with the equation of motions and ([551) we find the same 
expression for the energy density, p, and the pressure, p, with Eqs. (|26l) and 
(|27p. respectively. Thus the Friedman equations read 

3 - A = SirG ( n <#) + a„(V^> 2 " J , (42) 



n=l 
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and 



(*) 



a(t) 



At) 

8ttG ( ]T a„(2 



A 



(43) 



As is mentioned in the previous section, the Eq. f)43[) can be derived from the 
Eq. (|42l) with the help of continuous equation for the homogenous pressure and 
density from the cosmic fermionic components. 
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Figure 6: The behavior of the Hubble and the deceleration parameters for 
A/(87rGm 4 ) = 0.1, (iptp)(t = 0) = m 3 ,a{t = 0) 
/3 = 1, /?im 6 = 0, 1, 2 and /3 m = for m ^ 0, 1. 
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Figure 7: The behavior of the Hubble and the deceleration parameters for 
A/(87rGm 4 ) = 0.1, (inp){t = 0) = m 3 ,a(i = 0) = m" 1 , aim 2 = 10, n max = 1, 
/3 = 1, /?im 6 = -1, /3 2 m 12 = 1, 2, 3 and P m = for m ^ 0, 1, 2. 

We again assume that the non-trivial solution is realized for some of spinor 
fields and numerically solve the Friedman equation (1421) under the non-trivial 



11 



solution pt)l) . The results depend on the parameters a n and f3 n . Varieties of 
evolution can be reproduced by choosing these parameters. Here we set finite 
values for a\ and /3 n with n 6 {—1,0,1,2}. To avoid the upper limit for the 
scale factor we introduce a cosmological constant. 

As is seen in Figs. [6] and [7] the acceleration period, q < 0, appears later 
for a larger (3± and p2 ■ If we introduce a negative «2 with a positive a^, a 
similar behavior with Fig. @] is observed for a small a. Therefore an early-time 
acceleration can be induced by the negative 0:2 . Non-canonical kinetic terms can 
control when the late-time acceleration period starts with the model parameters 
tuning properly. 




am 



(a) H 2 /{%nGm A ) ( b ) 9 

Figure 8: The solid lines show the behavior of the Hubble and the deceleration 
parameters for A/(87rGm 4 ) = 0.1, (4>ip)(t = 0) = m 3 ,a(t = 0) = m _1 , aim 2 = 
10, n max =l,/3 = 1, fi-imr® = -0.05, -0.1, -0.15 and (i m = for m ^ -1, 0. 
The dashed lines represent the behavior for A = 0. 

The term, /3_i(^^) -2 , raises up the Hubble parameter. As is shown in 
Fig. [SI an open universe can be realized without any cosmological constant like 
gravitational source for a negative /3_i. 

4 Condensation of spinor fields 

In QCD like theory the gauge interaction becomes stronger at low energy. The 
strong gauge interaction induces the condensation of spinor fields below a critical 
scale, a cr or a critical temperature T cr , see for example Ref. |40) . After the 
condensation the expectation value for the composite operator constructed by 
fermion and anti-fermion, -fiif), develops a non-vanishing value. It is known as the 
spontaneously chiral symmetry breaking result in QCD |36) . Thus it is natural 
to use the trivial solution (ipiji) = before the universe arrives the critical scale, 
a(t) < a cr . 

In a static background metric the multi-fermion interaction is evaluated in 
Refs. (4TJ |42j [43] ■ In the FRW universe the critical scale have to be fixed 
by observing the effective action of the theoretical model. For this purpose the 
quantum field theory in an unstable background metric is necessary. It seems to 
have many difficulties 44, 45^. We take the remaining problem to fix the critical 
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scale as future works. However, the symmetry breaking changes the solution 
for the equation of motion. It has played a decisive role for the evolution of the 
universe. Here we set the critical scale by hand and discuss the contribution 
from the symmetry breaking effects at certain scales. 

We suppose that the first order phase transition takes place, the solution for 
the equation of motion suddenly jumps from the trivial to the non-trivial one. 
Under the trivial solution the Hubble parameter is constant, which corresponds 
to the well known de Sitter universe phase. A static universe realized without 
any cosmological constant like terms. The cosmological constant can exponen- 
tially expand the universe. After the universe arrives at the critical scale, the 
expectation value for the composite operator ^tjj decelerates the cosmic expan- 
sion. 
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Figure 9: Behavior of the Hubble and the deceleration parameters for a cr m = 
1.5, A/(87rGrn 4 ) = 0.1, C = —1, aim 2 = 5 and n max = 1. 

We set the critical scale at a cr m = 1.5 by hand and draw a typical behavior 
of the Hubble and the deceleration parameters for the same parameter setting 
with a line for aim 2 = 5 in Fig. [5J As is seen in Fig. [5J the Hubble and the 
deceleration parameters are constant, H 2 — A/3 and q = —1 below the critical 
scale a < a cr . We observe the same solution with Fig. [3] above the critical 
scale. In this case the spinor field contributes in the intermediate scale. The 
cosmological constant dominated universe can be realized at the early and the 
late stage of the universe evolution. We note that the cosmological parameters 
moves smoothly at the critical scale if the phase transition is of the second order. 

5 Conclusions 

We have investigated the self-interacting spinor fields as a gravitational source 
for space-time evolution via the Einstein's gravity. It has been assumed that the 
background metric is homogeneous and isotropic in the large scale. We have 
solved the Einstein-Dirac equations for the spinor bilinear, ipip, in the FRW 
metric. There are two solutions. One of the solutions, 'tjjip — 0, has played no 
significant role for the space-time evolution. 

First, we have evaluated the non-trivial solution and showed some behaviors 
of the cosmological parameters. It has been observed that a self-interacting 
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potential can induce the early time cosmic acceleration. It is consistent with re- 
sult in Ref. [3] . The expectation value for the mass term dominates the vacuum 
energy at a late stage of the universe evolution. Because of a negative expec- 
tation value for the spinor bilinear operator, tptp, mass term decelerates cosmic 
expansion. Only a closed universe can be realized without another gravitational 
source. In this work we have introduced a cosmological constant term. The 
cosmological constant can promote the accelerated expansion of the Universe at 
a late stage. It is shown that the non-canonical kinetic term controls when the 
late time acceleration starts. 

Non-vanishing expectation value for the spinor bilinear breaks the chiral 
symmetry for the spinor field. It is believed that high temperature, high density 
and strong curvature may restore the broken symmetry at the early universe. 
Thus it is not always valid to adopt the non-trivial solution of the Einstcin-Dirac 
equations. The trivial solution is realized before the critical scale, a < a cr . Here 
we have set a critical scale, a cr , by hand and demonstrate the possibility that 
the spinor fields can decelerate the cosmic expansion in an intermediate scale. 

It should be noted that the late time accelerated expansion or the dark energy 
phenomena can be produced without a cosmological constant. An interesting 
mechanisms is found in a modified gravity. For general review of the later time 
accelerated expansion in the modified gravity, see Ref. It is not hard to 

take account of spinor source in above scheme. 

In this paper we restrict the potential as a function of only the scalar in- 
variant operator, tpip, and consider isotropic universe evolution for large scale. 
Other types of spinor bilinear may generate an angular momentum. It is known 
that an anisotropic universe is generally realized due to the spinor angular mo- 
mentum. Hence, it is interesting to include other types of spinor bilinear and 
study the model in an anisotropic universe. We have not calculated the critical 
scale, a cr , starting from the action. It should be fixed by observing the sta- 
tionary point of the effective action. To evaluate the effective action we should 
take into account the quantum and thermal effects in an expanding universe 
background. We will continue our work further and hope to report on these 
problems. 
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